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Abstract. A suppression in the spectrum of ultrahigh-energy (UHE, > 10^* eV) 
neutrinos will be present in extra-dimensional scenarios, due to enhanced neutrino- 
antineutrino annihilation processes with the supernova relic neutrinos. In the n > 4 
scenario, being n the number of extra dimensions, neutrinos can not be responsible 
for the highest energy events observed in the UHE cosmic ray spectrum. A direct 
implication of these extra-dimensional interactions would be the absence of UHE 
neutrinos in ongoing and future neutrino telescopes. 
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1. Introduction 

Experimental high-energy neutrino astronomy is developing very rapidly. There 
exist a number of experiments (AMANDA II [1], RICE ^, ANITA [3], Icecube 
ANTARES [5]) that are currently analyzing or starting to take data. In the future there 
are planned projects (ARIANNA [6], AURA, NEMO, ACORNE) that will benefit from 
improved detection techniques and larger effective detection volumes. 

A guaranteed source of UHE neutrino fluxes are the so-called cosmogenic GZK 
neutrinos, which are originated by the interactions of extragalactic UHE cosmic ray 
(CR) protons with CMB photons dominantly via A"*" processes and subsequent charged 
pion decays. Cosmogenic neutrinos are typically characterized by a spectrum peaking 
in the 10^^~^^ eV energy range, depending on the redshift of the CR sources. Ongoing 
and future experiments expect to detect a few GZK neutrino events; the precise number 
depends on the full exposure of the instruments as well as on the production model. 
Direct emission of UHE neutrinos from the CR sources is expected but uncertain. Decays 
of topological defects or supermassive particles, leftover fossils from the GUT era, is 
speculative. Nevertheless, both mechanisms would produce neutrino fluxes with energies 
comparable to or higher than those associated to the GZK fluxes. These neutrinos could 
interact with 1.95 °K CMB neutrinos (CuB) via the standard model (SM) reaction 
up — > provided that they are extremely energetic (10^^~^^ eV) [3 El [9l [10]. We do 
not explore these speculative neutrino fluxes in the present study. 

In this study, we focus on the depletion of the GZK cosmogenic neutrino fluxes via 
strongly interacting annihilation processes with other neutrino relics that also permeate 
the universe: the diffuse supernova relic neutrinos (DSNz/), that represent the flux of 
neutrinos from all supernova explosions that occurred during the universe's history. 
The DSNi/ direct detection is still elusive. The most stringent experimental current 
limit to the DSN relic z/g flux is 1.2 cm~^s~^ at 90% CL, from the SuperKamiokande 
experiment |llj. The presence of strongly interacting processes, such as the exchange 
of massive spin-2 particles in theories of large extra-dimensions [T2 | ITS f IT4j. can modify 
the I'l' annihilation cross section. This effect would take place at high values of the 
squared center-of-mass energy s, yielding a annihilation cross section that is larger 
than the cross section for the SM process vv^^^ ^ . In principle, the UHE cosmogenic 
neutrinos can annihilate with both the Ci/B |T5] and DSNz/ via extra-dimensional 
enhanced cross sections, which we discuss next. 

2. Neutrino annihilation in extra-dimensional models 

We consider the following annihilation cross sections for n extra dimensions p3l [T5] 




(1) 
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respectively to produce KK gravitons, fermion- and 7- pairs. Here JF^ = vr^ + 
4:P{Ms/ ^/s) and we use I{Ms/\/s) as given in Ref. [H]. The "new physics" scale 
Ms is constrained from astrophysical considerations such as star cooling by graviton 
emission [16] and from collider searches ^17|. In particular, we use Ms = 701 TeV, 
25.5 TeV and 2.77 TeV for n = 2, 3 and 4, the most stringent current constraints from 
heating of neutron-stars [18]. For n > 4, the most stringent lower bounds are from the 
DO collider experiment at the Tevatron, which sets the 95% CL limits for n=5, 6 and 
7 equal to 0.97 TeV, 0.9 TeV and 0.85 TeV, respectively [I?]. In the n = 5 scenario, 
the total annihilation cross section is ~ 4 x 10^^^ cm^ at ^/s ^ 14 TeV, which 
roughly corresponds to a 10^^ eV GZK neutrino interacting with a 10 MeV DSN relic 
antineutrino. The cross section quoted above is therefore many orders of magnitude 
larger than the SM cross section o'^pL.au — 8 x 10~^^ cm^ at the same ^/s and scales as 
~ s^ for n = 5 and s ^ Ms- 

The neutrino interactions in Eqs. ([1]) are independent of the neutrino flavor. Brane- 
bulk couplings are flavor blind and consequently the exchange of the KK gravitons 
is unaffected by the electron, muon or tau nature of the DSN (anti)neutrinos, except 
corrections proportional to the squared mass sphttings divided by s, which are negligible 

(0(10-27)) g. 

A word of caution is needed here regarding the extra-dimensional scenario, which 
is an effective theory valid for s ~ M|. At some energy scale s ~ M|, this theory 
is supposed to match onto a more fundamental theory of quantum gravity. It is not 
known how to do this matching. A phenomenological approach is to assume that the 
neutrino interaction cross sections in the s ~ M| energy range behave similarly to the 
cross sections in the s ~ M| energy regime, up to some cutoff A. The value of A is 
presumably somewhere between Ms and -Emax? where the latter is the scale at which 
perturbative unitarity would be violated [13]. For the models we consider -Emax is always 
greater than 5.6M5. 

Within the context of extra-dimensional models, the vN cross sections will be 
enhanced as well [IHl EDI EH [22], [23l [211 EH] , providing a possible explanation for the 
events above the GZK cut-off as explored in Refs. [261 HSl EH 128]. However, as we will 
discuss shortly, 10^° eV neutrinos would annihilate with DSNi/ on their flight to the 
Earth rather than producing an extended air shower in the atmosphere, via enhanced 
vN cross section, in the large extra-dimensional models. The advantage of exploring 
the vv annihilation channel is that extradimensional signatures would occur at lower 
energy, compared to the signatures in the commonly explored vN interaction. 

X This flavor blindness character of the extra-dimensional model presented here no longer holds if one 
or more of the neutrino species are in the bulk. Such a possibility is not considered through the present 
discussion. 
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3. Supernova relic neutrino density and UHE neutrino propagation 

A number of authors have predicted the DSNi/ flux. For a recent appraisal of the 
theoretical and computational status, see Ref. [29] and references therein. Here we 
follow closely the derivation given in Ref. [30]. A fit to the neutrino spectra from 
numerical simulations of a SN is [31l |32] 



dE, T{l + (3,)El \EJ ' ^ > 

where the average energy E, = 15.4 MeV and 21.6 MeV respectively for z/g and 
corresponding to all other non-electron anti-neutrino and neutrino flavors. The 
spectral indices are = 3.8 and = 1.8 while the total neutrino energies are 
Lp^ ~ Lj,^ = 5 X 10^2 erg. For u^, we use E,^ = 11 MeV [323, ^i.^ - Lp, and f3,^ = f3p^. 
Neutrino conversion inside the star mixes the different neutrino flavors and therefore 
the relic (anti) neutrino flavor spectra at the stellar surface will differ from the original 
ones. The final flavor spectra will depend on the neutrino mass ordering (normal versus 
inverted) and the adiabaticity of the transitions in the resonance layers, see Ref. [33] for 
a complete description. As we will explain further below, the interactions we explore 
here are flavor blind and therefore the GZK (anti) neutrino will interact with the three 
(neutrino) antineutrino flavors. Therefore we do not need to account for conversion 
effects and the relevant quantity would be the total antineutrino (neutrino) SN relic 
neutrino spectra, given by: 

dNp^^^K^^o^ m 
dE, dE, dE,' ^ ^ 

that is, the sum of the three flavor spectra. 

The redshift-dependent SN rate is a fraction 0.0122Mq^ of the star formation rate 



and is given, e.g. SFl model in Ref. [34], by 

exp(3.42;) 



R,„(z) = 0.0122 X 0.32/1 



™exp(3.8z) +45 

with a Hubble constant Hq = 70hm km s~^ Mpc~^ and ACDM cosmology. The other 
parameters are Qm = 0.3 and = 0.7. The differential number density of SN relic 
neutrinos at present from all past SNe up to a maximum redshift -Zgn.max is then [30] 

dnpi^,) /•2sn,max dt .dNp(,^ 

Here {dt/dz)~^ = -Ho{l + 2)[^]^(l + zf + ^aI^/^ and E, = E'J{1 + z) is the redshift- 
corrected observed energy. 

While the number density of the DSNi/ (10^^ cm~'^ for the sum of the three 
(anti) neutrino flavors) is orders of magnitude smaller than those for the Cz/B relics 
(56 cm~^ per each (anti) neutrino flavor), the average energy of the DSNi/ is tens of MeV, 
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compared to the 10~^ eV for CuB relics. Therefore, the UHE neutrino mean-free-path, 
mfp = 1/ CTupi^u is many orders of magnitude smaller in the case of the less abundant, but 
more energetic DSNz/ compared to the Cz/B relics. If the strongly interacting processes 
deplete the UHE cosmogenic neutrino fluxes, the dominant attenuator will be the DSNi/ 
targets, which we discuss more quantitatively below. 

An UHE u of observed energy -Ej^^uhe niay interact with a DSNi/ at redshift z' on 
its way via processes in Eq. ([1]) and annihilate. The corresponding s ^ 2£',^uhe(l + 
z')Ey^s^{l + z), ignoring the u masses. We use the maximum SN u energy to be 
E'„ rr,^^ = 60 MeV in the SN rest frame. The inverse mfp for uu annihilation is then 



-sn.max 



fit 

C-\E,^^^,-z')= I ' dz-{l + z)Rs 



dz 



X 



i^.sn.max 



■■ ' ■ - sn 



dEl^sn-T-^(^uu{s). (6) 

The mfp for a 10^^ eV neutrino to annihilate with a DSNz/ via the SM process v^^^n 
is 10^^ Mpc, which exceeds the Hubble distance. Within the n = 5 extra-dimensional 
model, the annihilation cross section is greatly enhanced at high energies, and the mfp 
for a 10^^ eV neutrino is ~ 12 Mpc in our local universe {z' ~ 0), which is less than 
the GZK radius. Even for the = 4 extra-dimensional model, the mfp for the highest 
energy CR, 3 x 10^'' eV, is ~ 127 Mpc which is comparable to the GZK radius. To 
explain GZK CR data with UHE neutrinos through enhanced vN cross section requires 
> 4. Thus UHE neutrinos propagating from outside the GZK radius can not be the 
candidates for GZK CR events, since they would be absorbed by DSNz/. 

We can now calculate the survival probability for an UHE z/ created at redshift z^he 
to reach Earth as 

rZuhe ^ dt 



PiE^,uhc; Zuhc) = exp 



-c 



/ dz'-^C-\E,,^^,;z') 
Jo dz' 



exp 



C f^nhc dz' 

-K, 



X 



-sn.max 



^0 ^0 (1 + z') yru(rT^TTa^ 

dz exp(3.4z) 



1 + 2)3/2 exp(3.82) + 45 



X / dE^^sn-T^ 

W-CV;/ sn 



(7) 



where ICc/Hl ^ 2.45 x lQ~^^hj^ cm and the differential SN z/ spectrum is 
dNp^s^/ dE^s^ ^ lO''^ MeV~^. Large z/z/ cross section then suppresses UHE neutrinos. 
We discuss UHE z/ fluxes that will be attenuated by z/z/ annihilation next. 

4. Ultrahigh-energy neutrino flux 

The CR energy generation rate per unit volume in our local universe in the energy range 
]^gi9-2i gY jg 5 X lO'^^ erg Mpc~3 yr"^ [35]. Assuming an injection spectrum for 



Ultrahigh- energy neutrino flux as a probe of large extra- dimensions 



6 



CR protons dNp/dE^ oc E^^, as typically expected, we define a convenient conversion 
formula 

■'VCR 



47ri7oln(102VlOi9) 
^ 7.1 X lO-^/ifo^ GeV cm-2 sr^S (8) 

which is proportional to the CR flux E^Jp above 10^^ eV. We will use Eq. (jH]) to fix 
the normalization of UHE v fluxes. The CR sources may also evolve with redshift as 
S{z) = (l + zfhiz < 1.9, (1 + 1.9)3 for 1.9 < ^ < 2.7 and exp[(2.7-2)/2.7] for 2 > 2.7 
[35]. 

The Waxman-Bahcall (WB) bound on UHE u flux [35] is based on CRs that interact 
at their sources and lose all their energy equally to charged and neutral pions. The 
resulting z/^ flux is given by 



/fi^(l + Zuhc)^ + 



after integrating over CR source evolution and uu annihilation probability in Eq. ([7]). 

If UHE CRs interact with CMB photons in the local universe then the resulting 
GZK neutrino flux would be 

E^Mz ~ 0) cx AfcK J dEl'^Y{El, E,, z r-. 0) (10) 

Here Y is called the neutrino yield function as in Ref. [SB] and is the number of secondary 
neutrinos generated per unit energy interval by a CR proton of energy E^. We use a fit 
to Y{Ep, Ey, 2; ~ 0) corresponding to z/^ and z/^ from a CR proton propagating 200 Mpc 
as generated by the SOPHIA Monte Carlo code as reported in Ref. [36j. The GZK 1/ 
spectra are fully evolved by 200 Mpc in our local universe and over smaller distance at 
higher redshift. Our calculation shows that this distance is much shorter than the mfp 
for vN interactions of UHE CRs with DSNz/ in n > 4 large extra-dimensional models. 
Thus we calculate the effect of vv annihilation assuming that a fully evolved GZK u 
flux exist at a given redshift of interaction. 

The GZK v flux integrated over all CR sources, after taking into account the redshift 
evolution of the neutrino yield function Y{E^, E^, z) = Y{Ep{l -\- z),Ei,{l + z)^, z ~ 0) 
[36] . the source evolution S{z) and finally the survival probability P{Ey; z^he) in Eq. ([7]), 
is given by 

S{z^^,)P{E 



EyJu,GZK = A/cR / C^^uhc 

Jo 



fim(l + 2;uhc)3 + 



X I dE;^YiE;,E,,z^^,). (11) 



In case of no vP annihilation, P{Ey] Zuhe) = 1 and the flux is the same as in Ref. 

We have numerically evaluated the GZK flux, both without and with 
annihilation, using ^max = ^uhe = 2;sn,max = 5 and in the energy range 10^^ eV < E'^ < 
10^^ eV with an exponential cutoff of the cx E~'^ spectrum at 3 x 10^^ eV as in Ref. |36j . 
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Figure 1. UHE h' fluxes from the cosmic ray protons interacting at the source (WB) 
and in CMB (GZK). ff vD annihilation is important, as in large extra-dimensional 
models (shown here for n = 5 case with dotted curves labeled 1 and 2), then UHE v 
fluxes would be suppressed. Also shown are the projected sensitivities for the ANITA 
(50 days) and the proposed ARIANNA (6 months) UHE neutrino experiments at the 
South Pole. 

The results for the GKZ cosmogenic z/^ flux are depicted in Fig. [T], assuming a n = 5 
extra-dimensional scenario (the dotted curves labeled 1 and 2 corresponds to A = oo 
and 20 TeV respectively; allowing the cross sections in Eq. ([T]) to grow below y/s = A 
and become fiat above). Also shown is the WB flux without and with uu annihilation. 
Notice that the n = 5 extra-dimensional scenario leaves a clear imprint on the GZK 
cosmogenic neutrino fluxes, which would be abruptly truncated above E > 10^'' eV. 
This characteristic feature in the GZK cosmogenic fluxes could be recognized by the 
presence of a dip in the neutrino spectra, provided the detection technique has a low 
enough energy threshold. For ongoing and future UHE neutrino experiments with higher 
energy thresholds {E > 10^^ eV), such as ANITA and ARIANNA shown in Fig. [H 
there would be an absence of neutrino induced events caused by strongly interacting, 
KK-modes mediated processes. For the n < 5 extradimensional models, the UHE 
neutrino flux suppression would occur at UHE neutrino energies E > 10^^"^° eV, where 
the cosmogenic neutrino fluxes are smaller and consequently, also the statistics expected 
in ongoing and future UHE neutrino observatories would be reduced. 

Figure [2] depicts the the GZK cosmogenic flux with and without extradimensional 
suppression for the case nature has ra = 4, 5, 6 and 7 extra dimensions. For n < 4, the 
UHE neutrino flux suppression is subtle and therefore it would be highly challenging 
and difficult to detect experimentally. Figure [3] illustrates the WB flux without and 
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Figure 2. The solid line depicts the UHE v fluxes from the cosmic ray protons 
interacting with CMB photons (GZK neutrino fluxes). The dotted, long-dashed, dot- 
dashed and short-dashed curves illustrate the GZK neutrino fluxes for the case of 
n = 4, 5, 6 and 7 extradimensions exist in nature, respectively. 

with vv annihilation. If n < 5, tracking the extra-dimensional induced suppression dip 
would be more difficult in general. Note that an increase of vN cross section, expected 
in this scenario, do not significantly increase the detector sensitivity because of a steeply 
falling z/ flux and a decreasing angular acceptance with increasing energy (see, e.g., [2]). 

5. Summary and conclusions 

We have shown that UHE neutrinos will be absorbed, in theoretical models that predict 
fast-rising cross sections such as large extra-dimensional models, by a diffuse background 
of 10 MeV neutrinos provided by all core-collapse SNe in the history of the universe. 
Detection of neutrinos from the SN 1987A proves the existence of such neutrinos, and 
upcoming megaton detectors will measure the diffuse flux to a good accuracy. 

If there exist n > 5 large extra-dimensions in nature, and the DSNz/ flux is detected 
at the level of the current theoretical models, then UHE neutrinos can not be the 
primaries of the super GZK events, since the UHE neutrino fluxes will suffer a cutoff^ in 
their energy spectra in the 10^^~^^ eV energy range. On the other hand, a detection of 
GZK neutrinos at energies E > 10^^ eV could imply the absence of n > 5 large extra- 
dimensions in nature, and therefore eliminating such models. For tt, < 5 extradimensions, 
neutrinos could be the UHE CR primaries if the vN cross-section is sufficiently enhanced 
to mimic hadronic cross- sect ion. 

In case the DSNz/ flux is detected at a much lower level, then the dip in the UHE 
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Figure 3. Same as Fig.[2]but for UHE ly fluxes from the cosmic ray protons interacting 
at the source (WB neutrino fluxes). 



neutrino spectrum, due to absorption by DSNz/, would be shifted to higher energy. Note 
that uu annihilation by UHE neutrinos would not produce 7-rays over the EGRET limit, 
since the primary UHE CR interactions with CMB and infrared photons can not account 
for the observed diffuse 7-ray flux [37] . Also the GZK CRs are not affected due to large 
uN cross section, since they are expected to be produced within ~ 50 Mpc, a radius 
smaller than the uN mfp with enhanced cross section. 

Measuring an enhancement of UHE neutrino cross sections at ongoing or future 
neutrino observatories, will be therefore extremely difficult, since in these scenarios 
the GZK cosmogenic neutrino fluxes would be depleted in their way to the Earth via 
annihilation with the DSNi/ background. 
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